Polarization properties of excitonic qu-bits in single self-assembled quantum dots 
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(Dated:) 

We investigate polarization properties of neutral exciton emission in single self-assembled InAs/GaAs quan- 
tum dots. The in-plane shape and strain anisotropy strongly couple the heavy and light hole states and lead 
to large optical anisotropy with non-orthogonal linearly polarized states misaligned with respect to the crys- 
tallographic axes. Owing to a waveguiding experimental configuration, luminescence polarization along the 
growth axis has been observed revealing the presence of shear components of the deformation tensor out of the 
growth plane. Resonant luminescence experiments allowed determining the oscillator strength ratio of the two 
exciton eigenstates. Valence band mixing governs this ratio and can be very different from dot to dot, however 
the polarization anisotropy axis is quite fixed inside a scanned area of one jim 1 and indicates that the in-plane 
deformation direction to which it is related has a correlation length of the order of magnitude of one jim 1 . 

PACS numbers: 71.35.-y, 78.47.Fg, 78.67.Hc 



I. INTRODUCTION 

Addressing and manipulating an individual spin in 
a single quantum dot (QD) is a fundamental step to 
achieve for eventual information processing based on 
quantum mechanics. 1-3 Neutral excitons can also be 
used as single qu-bits but their polarization proper- 
ties have to be precisely known in order to address a 
specific state. It is well established now that for self- 
assembled quantum dots, in-plane shape anisotropy in- 
duced during growth leads to two linearly polarized 
orthogonal exciton states, splitted by the exciton fine 
structure splitting (FSS), 4-11 having in principle the 
same oscillator strength. 

However, recent experimental observations show 
that the linear states are neither aligned along the crys- 
tallographic axes, nor orthogonal between them and 
an important linear polarization degree can be found. 6 
These facts strongly indicate that valence-band mix- 
ing (VBM) due to the anisotropic confinement poten- 
tial and strain effects should be taken into account to 
explain the strong polarization anisotropy observed in 
the emission of individual QDs, as well as the exciton 
eigenstates orientation. 

Indeed, in II-VI CdTe/ZnTe 1213 or III-V 
InAlAs/AlGaAs 14 QDs, VBM has been related to 
the strain anisotropy induced during growth pro- 
cess. In strain-free GaAs QDs grown by droplet 
epitaxy it has been shown that VBM exists due to 
the anisotropy of the confinement potential, 15 while 
in strained self-assembled InAs/GaAs QDs, VBM 
effects have been investigated in order to estimate the 
exciton spin relaxation time, 16 its influence on the 
exciton FSS, 17 ' 18 or the hyperfine coupling. 19 To our 
knowledge, no experimental study includes both strain 
effects and confinement anisotropy in order to explain 



the observed results, which can be fully described 
by combining the Luttinger-Kohn and the Bir-Pikus 
Hamiltonians. 20 ' 21 

We report on polarization resolved micro- 
photoluminescence (jiPL) of the neutral exciton 
emission in single self-assembled InAs/GaAs QDs. 
We investigate the role of VBM on polarization 
properties, due to confinement potential anisotropy 
and local strain effects described by a Luttinger-Kohn 
and Bir-Pikus Hamiltonian. 

Since strain effects and dot shape anisotropy depend 
strongly on growth conditions they vary from sample to 
sample and from dot to dot. Therefore, a large disper- 
sion of values for the heavy to light-hole mixing param- 
eters can be found in the literature. 13-15 ' 18 Moreover, 
since morphological details and precise values for the 
different strain components can not be obtained at the 
same time for each studied QD by optical experiments, 
it is quite impossible to know precisely the origin of 
valence-band mixing. 

Nevertheless, from the study of polarization 
anisotropy on a very large statistical ensemble of indi- 
vidual QDs (about 300 dots on two different samples), 
we can conclude that: (i) the polarization anisotropy 
axis direction is dominated by strain-induced VBM, (ii) 
the polarization main axis is not aligned along the crys- 
tallographic axes but has a given direction within one 
probed micron, and (iii) the VBM strength can be large 
and can vary significantly from dot to dot in the same 
one probed micron. These results seem to indicate that 
if the polarization main axis is related to the in-plane 
main strain orientation, 13,14 the mean components of 
the strain tensor are quite homogeneous for all dots in- 
side a Ijlm 2 area. However the strength of the mixing 
seems to be mostly related to confinement anisotropy. 

Using an original experimental configuration where 
the QDs layer is embedded in a one-dimensional 



waveguide, resonant emission of the neutral exciton in 
single QDs could be analyzed and yield the two eigen- 
states oscillator strength ratio. Using this geometry, 
polarization-resolved jlPL analyzed along the growth 
axis could also be recorded, revealing coupling be- 
tween same-spin heavy and light hole states due to off- 
plane strain components. In conventional geometries 
this can not be evidenced and only in-plane strain com- 
ponents are usually taken into account. 13 ' 14 

The paper is organized as follows: in section II we 
give some details on the samples and the experimental 
configuration. Section III is devoted to the in-plane po- 
larization - resolved jlPL experiments on single QDs 
excited non-resonantly (part A) and resonantly (part 
B). In parallel, simulations have been performed in or- 
der to give the general trends of the polarization prop- 
erties influenced by mixing effects. In order to deduce 
the amplitude of the hh-lh coupling and the anisotropy 
axis of the polarization the experimental results are fit- 
ted taking into account the VBM. Moreover, an anal- 
ysis of the eigenstates polarization misalignment with 
respect to the crystallographic axes is made. Section IV 
presents the results of polarization resolved jlPL along 
the z-axis which allowed to evidence the shear strain 
components out of the growth plane. 



II. SAMPLES AND EXPERIMENTAL SET-UP 

InAs/GaAs self-assembled QDs, were grown by 
MBE on a planar [001] GaAs substrate. In order to 
address the fundamental exciton state at resonance, 
the dots are embedded in a one-dimensional GaAlAs 
monomode waveguide (WG). 22 The wave-guiding ge- 
ometry has several advantages: first, due to the spatial 
confinement of light, the volume of the optical mode 
is reduced and the coupling with the QDs is enhanced. 
Second, contrary to the microcavity case, there is no 
need to match the cavity mode energy with the QD 
emission energy in order to achieve a strong coupling 
between light and matter. Finally, the laser light prop- 
agation axis can be perpendicular to the luminescence 
detection direction. Indeed, while the laser propagates 
in the (x,y) dots plane, the single QD luminescence 
can be collected from the WG top surface by a con- 
focal jlPL detection set-up (see Fig. 1). Since the 
laser beam is confined in the guided mode, the scat- 
tered light is greatly suppressed and the resonant lu- 
minescence is almost background-free, at low pump 
power. Using this geometrical configuration, resonant 
Rabi oscillations of the exciton state in a single QD 
have been observed. 23 Coherent control experiments 
have also been performed in order to determine the 
main decoherence mechanisms in the system. 24 More- 
over, in this geometry, it is also possible to excite a dot 
by the sample top surface and to detect the PL at the 
edge of the WG. This last configuration allows us to 
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Figure 1: (Color on-line) Waveguiding configuration. The 
laser propagates along the x-axis and is polarized along y. 
The jJ.PL is detected along the z-axis and the polarization is 
analyzed by a half- wave plate followed by a polarizer. The 
WL and the QDs are represented as a red plane embedded 
in a waveguide formed by AlGaAs/ GaAs layers in blue and 
grey respectively. One micron ridges are etched on the sam- 
ple surface to form the one-dimensional monomode optical 
guide. 



investigate the PL properties along the z-axis, which in 
usual back scattering geometry can not be resolved. 

The QDs are excited by ps pulses provided by a 
mode-locked Ti: sapphire laser and focused on the sam- 
ple using a microscope objective. The PL signal is 
then collected by another objective (N.A. = 0.5) and 
detected with a 0.5 m grating spectrometer and a liquid- 
nitrogen-cooled charge-coupled device camera. All the 
objectives are mounted on piezo-electric actuators so 
they can easily control the coupling to the WG and 
move along the ridge to probe different QDs. 

The spatial resolution of the set-up is diffraction- 
limited to 1/lm 2 , while its spectral resolution is about 
A5jleV. In order to resolve the neutral exciton FSS, a 
Fabry-Perot interferometer, used as interferential filter 
can be inserted in the collection path, improving our 
spectral resolution to 0.6/1 eV. Then we were able to 
measure the exciton FSS Spss* which varies depending 
on the dot from 20 to 145 jleV . The laser is polarized 
along [110] and the sample is fixed on the cold finger 
of a three-axes optical cryostat (10K) specifically de- 
signed for the wave-guiding geometry. 

Two experimental configurations can be used for 
the polarization-resolved experiments: either the pump 
laser is guided and the jlPL is collected from the top 
surface of the WG, or the dot is excited from the top 
surface and its PL is collected at the edge of the WG. 
In both cases, the PL polarization is analyzed using a 
rotating half-wave plate (HWP) and a fixed linear po- 
larizer, in order to avoid selection from the polarization 
response of the whole detection setup. 
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III. VALENCE-BAND MIXING EFFECTS IN THE 
(X,Y) DOTS PLANE 

In the case of QDs with revolution symmetry around 
the growth direction z ([001] axis), the two neutral exci- 
ton bright eigenstates | + 1 ) and | — 1 ) , with a total angu- 
lar momentum / = ±1, are degenerated in energy and 
their emission is circularly polarized, left (<7 + ) or right 
(<7~) in the dots plane (x,y). When the system symme- 
try is reduced, the electron-hole exchange interaction 
lifts this degeneracy and the two bright eigenstates are 
the two linear combinations 7 ' 8 |X) = ^4= ( | +1 ) — | — 1 ) ) 

and |y) = 4j(|+l) + | — 1)), whose emission is linearly 
polarized along the crystallographic axes [110] (x-axis) 
and [110] (y-axis). 4-9 These two eigenstates are split 
in energy by the FSS (Spss = \Ex — Ey\), and under 
non-resonant excitation the emission intensity of both 
states is the same, whatever the difference in oscillator 
strength. 6 

If we now consider valence-band mixing (VBM) 
between heavy hole (HH) and light hole (LH) states 
which is due in general to in-plane shape asymmetry 
and anisotropic strain effects, the neutral exciton eigen- 
states become linear combinations of the elliptically- 
polarized bright exciton states |±1), which can be writ- 
ten as: 

|±l) = 0^S2| T i;±^)+/3 e ±2l >| T i; T i) (1) 

The first spin state corresponds to the electron spin 
projection, while the second corresponds to the heavy 
(±3/2) and light (±1/2) hole spin projection. The 
simulation adjustable parameters j3 and \j/ represent 
the amplitude and phase of the mixing as defined in 
the Luttinger-Kohn and Bir-Pikus (LKBP) Hamilto- 
nian (see appendix B). Emission of these two states 
can then be elliptically polarized and their main axis 
can be tilted by an angle \\f with respect to the y-axis 
[110]. The two neutral exciton eigenstates are then: 

\X) = -±=(\+l)-e 2id \-l)) (2) 

|?) = -L(|+l)+ e ^|-l)) (3) 

An additional angle 6 (see Fig. 2) has been introduced 
to take into account the QD main elongation axis ori- 
entation with respect to [1 10]. As discussed in Refs. 13 
and 14, this angle can be different from the crystallo- 
graphic axes and \jf. 

Depending on the LH ratio j3 and the angle i/a, the 
two eigenstates \X) and \Y) can be tilted with respect to 
the crystallographic axes [110] and [110], without be- 
ing perpendicular anymore. The tilt angles are defined 
as (f>x and §f respectively (see Fig 2(c)) and the angle 
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(a) 




(b) (c) 



Figure 2: Definition of the different angles involved in the 
simulations, (a) 6 is the main dot orientation axis, (b) \j/ is 
the main emission anisotropy axis and a is the angle between 
the analyzer axis and [110]. (c) (j) x and (j)f are the eigenstates 
tilt angles and Q. is the angle between the two \X) and \Y) 
eigenstates. 

between them is Q. In addition, under non-resonant 
excitation, because of the VBM, the two states have no 
longer the same emission intensity, resulting in the blue 
and red curves represented in Fig 2(c). 

A. Non-resonant excitation 

By calculating the dipole matrix elements for the 
transitions between the conduction-band and the Ts 
valence-band states (see appendix C), the two eigen- 
states emission intensity, passing through a linear po- 
larizer with an angle a with respect to the x-axis [110] 
(see Fig. 2(b)) can be determined. Under non-resonant 
excitation, assuming that the two levels are equally 
populated, independently of the laser polarization and 
their respective oscillator strength, we obtain the nor- 
malized intensity: 

Ix(a) = [ v / l-j3 2 cos(a±e) 

--^=cos(a-0±2y/ r )] 2 (4) 
v3 

I f (a) = [ v / l-j3 2 sin(a±e) 

±-^sin(a-e±2y/)] 2 (5) 




4 



= 0° 




(a) /3=0 



(b) p =0.25; v^ = 0° 




(c) p = 0.25; y/ = 45° (d) j8 = 0.25; y = 90° 



= +lO° 




(e) j8=0 



(f) P =0.25; v = 0° 



Figure 3: Calculated polar diagrams for |X) (blue curves) 
and |F) (red curves) emission intensity. The purple curve is 
the total emission intensity. The different parameters are: for 
6 = 0°, (a) p = (no VBM); (b) p = 0.25 and yr = 0°; (c) 
p = 0.25 and yr = 45°; (d) p = 0.25 and yf = 90°, and for 
6 = +10°, (e) P = (no VBM); (f) J3 = 0.25 and y = 0°. 



The total intensity in the (x,y) plan, which off- 
resonance is simply the sum of the two expressions 
above, is then written: 



4 r (a) = [l--/3 2 



-2J5^/^— ^cos(2(a- V ))] 



(6) 



We performed simulations using equations (4), (5) 
and (6) in order to understand the influence of the mix- 
ing parameters on the polarization properties of the ex- 
citon eigenstates. By plotting the total emission polar 
diagram (purple curve in Fig. 2(b)), the mixing param- 



eters j3 and yf can be determined. The \X) and \Y) polar 
diagrams give also the QD main axis orientation angle 
6 (see Fig. 3). Indeed, as shown in (6), the total inten- 
sity does not depend on 6 and the parameter yf is given 
by the polar diagram orientation angle with respect to 
[110]. /3 is given by the contrast of the total PL polar 
diagram. The \X) and \Y) polar diagrams can then be 
fitted to determine the unknown parameter 6. 

Calculated polar diagrams in Fig. 3 show that when 
VBM occurs (j3 ^ 0), the dot total emission is partially 
linearly polarized along a direction defined by yf, with 
a linear polarization degree (I Max - I min ) / \I Max + L) . 
In fact, this parameter depends only on j3 and can be 
expressed as: 



C(j3): 



2/V3(l-/32) 
3-2J3 2 



(7) 



Introducing, like in Refs. 13. and 18. the reduced 



variable 7] 



the linear polarization degree 
j|^2. Depending on/3 



takes the simpler form: C(tj) : 

and yf, the polarization of \X) (blue curve) and \Y) (red 
curve) emission can be tilted with respect to the crys- 
tallographic axes by an angle fa and fa respectively. 
Moreover, the VBM allows for the \X) and \Y) states 
not to be orthogonal. This is shown for instance in Fig. 
3(c), 3(d) and 3(f). The tilt angles fa and fa can be 
expressed as a function of tj, 6 and y/ by: 

fl -sin0 + T]sin(2vA-O) 



tan(^) 



COS0-T]COs(2l/A-0) 



sin0 + r]sin(2vA-0) 
cos0 + 7]cos(2i/a- 0) 



(9) 



First, let us consider the case where 6 = 0° (see Fig. 
3(a), 3(b), 3(c) and 3(d)), which occurs most often. The 
eigenstates tilting angles fa and fa only depend on the 
mixing parameters \j/ and j3 but in a non-trivial way. 
The general tendency is that for a fixed value of ft, fa 
increases, reaches a maximum value, then decreases as 
a function of y/. The larger j3 is, the larger is the maxi- 
mum value of fa. fa varies in an opposite way. 

Now, let us consider the case 0^0° (Fig. 3(e) and 
3(f)). For a given yf the presence of a small tilting of 
the dot changes the eigenstates orientation and thus the 
angle £1 between them. For instance, we can see that 
for y/ = 0\ the two eigenstates are perpendicular when 
9 = 0° (£} = 90°, Fig. 3(b)), but when = 10° the angle 
is £1 = 96° (Fig. 3(f)). Thus, it appears that £1 does not 
depend only on yf and /3 but depends also strongly on 

e. 

Polarization resolved jlPL measurements were per- 
formed on a very large statistical ensemble of indi- 
vidual QDs. Non resonant excitation energy has been 
used, typically in the absorption region of the wetting 
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(a) j8 = 0.15; y = 10°; = 0°; = 2°; 




(b) P = 0.15; y = 10°; = -6°; <fe = 9°; 
? = 3.5° 



Figure 4: Polar diagrams for |X) (bue dots) and \Y) (red dots) 
emission intensity, fitted by equations (4) (blue solid curve) 
and (5) (red solid curve) respectively. The total emission in- 
tensity is represented in purple. The two polar diagrams cor- 
respond to two different QDs located in the same detection 
spot with the same VBM parameters j8 =0.15 and if/ — 10°, 
(a) 6 = 0°, A = 954.5™ and S FSS = lOfdeV; (b) 6 = -6°, 
A = 929. %nm and 8 FSS = 125 neV. 



layer around 1.4 eV. Fig. 4(a) and 4(b) show two typ- 
ical experimental polar diagrams of two different QDs 
inside the same one jlm 2 detection spot. For these two 
particular dots, the same values for /3 and \j/ parameters 
were found (j3 = 0.15 and \j/ = 10°) but the orientation 
angle 6 was different (6 = 0° and 6 = — 6° for (a) and 
(b) respectively). Therefore, for a given total emission 
intensity (purple curve) it is not possible to determine 
the polarization direction of the two eigenstates (red 
and blue curves). In a general manner, we observed 
that all the dots (about 10 to 20) inside one flm 2 detec- 
tion spot have the same value of angle y/ but the VBM 
amplitude ft is very different for each dot, ranging be- 
tween and 0.5. In average, the mixing parameter is 
equal to 0.15 but with a large dispersion of ±0.11 (see 
Fig. 5(a)). The \j/ angle can also vary significantly de- 
pending on the dot position along the WG axis. How- 
ever the \j/ values are mainly concentrated around 10° 
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Figure 5: Histograms of (a) /3 and (b) y/ in one of the studied 
samples. The emission wavelength of the QDs range between 
900 nm and 940 nm. 



and 30° (Fig. 5(b)). Regarding the dots orientation axis 
0, rather small angles have been found, less than 15° 
and in most of the cases 6 was smaller than 5°. 

The fact that the angle \j/ is quite constant over a 
scale of one micron seems to indicate that y is related 
mostly to the in-plane strain anisotropy which is the 
only parameter for which we can define a correlation 
length. As also described in Ref. 13. 1// is correlated to 
the main axis of the in-plane strain field, suggesting an 
homogeneous strain field inside the detection spot. The 
important variations found for /3 on this same length 
scale indicate another physical origin having a strong 
influence on the mixing. Indeed, the VBM amplitude 
is related both to the strain and confinement potential 
anisotropy in the case of SK QDs. Since the in-plane 
strain anisotropy seems to be rather uniform over this 
distance, then the confinement potential should govern 
the amplitude of VBM. The non-uniformity of chemi- 
cal composition, 29 ' 30 the differences in shape and size 
of the dots 31 can create important fluctuations of the 
confinement potential resulting in very dispersed val- 
ues of the mixing parameter j3 . 

The strong VBM combined to the small misorien- 
tation of the dots results in tilted eigenstates with re- 
spect to the crystallographic axes. Tilt angles (j>x and 
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(a) P = 0.3; y = 25°; = 0°;/? = 0.7 




(b) = 3 (c) = 0.2 



Figure 6: Polar diagrams of the total emission intensity, for 
excitation off-resonance (purple curves) and on-resonance 
(pink curves), (a) represents experimental data for a given 
QD with p = 0.3, y = 25°, 6 = 0° and R = 0.7. (b) and (c) 
are calculated plots with the same parameters j8, \ff, 0, and for 
an oscillator strength ratio R = 3 and R = 0.2 respectively. 



§y reaching 20°, as well as an angle £1 ranging between 
75° and 90° were found. 

It is also worth noticing that in the case of the two 
dots presented in Fig. 4(a) and 4(b), the j3 parameter 
is the same but the emission wavelength and the FSS 
are different. Indeed, in a general manner we did not 
find any kind of correlation between the VBM ampli- 
tude and the emission wavelength or between the FSS 
and VBM, contrary to what has been recently reported 
in Ref. 18. This disagreement might be due to the dif- 
ferent confinement potential between the two kinds of 
samples, the emission energy of the dots being lower 
in energy in our case, suggesting a larger size. It is not 
also clear if the correlations reported in Ref. 1 8 are ob- 
served on a large area or in a very precise region on the 
sample. Our experimental results were obtained over 
hundreds of dots on different areas on the sample sur- 
face. 



B. Resonant excitation 

In-plane PL polarization under resonant excitation 
has also been investigated to determine the oscillator 
strengths ratio R = ff/fx between the two \Y) and 
\X) eigenstates. Although the excitation laser is lin- 
early polarized along the y-direction, both eigenstates 
can be excited since they can be tilted with respect to 
the crystallographic axes. Moreover the ps laser pulses 
are spectrally broader (« 500 jieV) than the FSS be- 
tween the two eigenstates (< 150/ieV). The pump 
power is kept low enough so the system is in the lin- 
ear regime. Taking into account different oscillator 
strengths for the two transitions, and the proper laser 
polarization, the neutral exciton emission intensity is 
always linearly-polarized . The total resonant normal- 
ized intensity reads now: 

I r (oc) = {[yl — j3 2 (cos0sin0£ -f sinflcos^y) 

— -^=(cos(0 — 2\jf) sin<^ -f sin(0 — 2\jf) cos (j)y)} 

(cos (j>x—R sin (j) Y ) cos a 

+[yl — j3 2 (— sin sin 0^ -f cosflcos^y) 

--^=(sin(0 —2\jf) sin<^ -cos(0 — 2\p) cos^y)] 
v3 

(sin0£ +7?cos0y) sin a} 2 (10) 

<j>x and (f)f can be replaced by expressions (8) and (9), 
so I r (oc) depends only on j3, \jf, 6 and R. Using the 
VBM parameters ft and \\f as well as the dot orienta- 
tion angle 6 that have been determined independently 
in non-resonant experiments, it is then possible for a 
given dot to obtain the ratio R. Fig. 6(a) shows a typ- 
ical polar diagram of the total emission intensity of 
a QD. The parameters /} = 0.3, \\f = 25° and 0=0° 
were found in off -resonance PL experiments and the 
ratio R = 0.7 was given by the orientation of the res- 
onant total intensity polar diagram (pink curve). Fig. 
6(b) and 6(c) show how the value of the ratio R can 
change the total intensity orientation although the off- 
resonance parameters are the same. If the eigenstates 
were aligned with respect to the crystallographic axes 
(\jf = 0° or 90° and 6 = 0° or 90°), then only the Y state 
would be excited. 

The oscillator strengths ratio was found to range be- 
tween 0.2 (fx > fy) and 1.5 (fx < f?), but we no- 
ticed that in most cases, fx was larger than fy. This 
ratio can be significantly different from the ratio be- 
tween the two eigenstates emission intensity found in 
off-resonance measurements. This fact confirms that in 
non resonant experiments the emission intensity is not 
proportional to the oscillator strength. 6 

Finally, knowing the eigenstates polarization and os- 
cillator strengths ratio for a given QD will allow ad- 
dressing a particular state with the proper polarization 
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and having a large oscillator strength. This will en- 
hance the possibilities of coherent manipulations of an 
excitonic qu-bit. 

The in-plane PL polar diagrams under resonant exci- 
tation can also allow to distinguish between the neutral 
exciton (linearly polarized) and the charged exciton (el- 
liptically polarized). This can be useful in non-resonant 
experiments, because if the FSS can not be resolved 
both neutral and charged exciton lines appear to be el- 
liptically polarized. 

IV. VALENCE-BAND MIXING ALONG THE 
GROWTH Z-AXIS 



I05« 




(a) (b) 



The original geometry of our set-up allows us, by ex- 
citing from the top surface and detecting the PL from 
the edge of the WG, to study the polarization prop- 
erties along the z-axis, which is impossible in usual 
backscattering excitation-detection geometry. 32 Obser- 
vation of PL with a polarization component along the 
z-axis in few dots, has led us to consider a mixing be- 
tween heavy and light hole states with the same spins 
due to the different off-plane shear strain components 
which are usually neglected in the LKBP Hamiltonian. 
The eigenstates are written as linear combinations of 
the elliptically -polarized bright exciton states: 

|±1) = A /l-jS2_ 7 2| T l ;± | ) 

+ e ±2 '>| T i ;T i) + re ±2 '%i;±i) dD 

where / represents opposite spin (in comparison to j8) 
LH amplitude coupling with t; the associated phase. 
The existence of 7 has no significant influence on the 
polarization properties for in-plane PL of dots, and 
does not influence the determination of the three other 
parameters j3, \jf and 6. Indeed, simulations (see Fig. 
7(a)) show that the presence of 7 changes neither the 
total emission intensity contrast (thus not changing j3 
and i/a), nor the two eigenstates emission intensity ori- 
entation, therefore not changing G. Only the normal- 
ization constants are slightly modified. We can also 
compare in Fig. 7(b) the polarization rates C(/3,7) = 

2 ^ w?fi2^ ~r^ as a function of the mixing parameter 
/3, for 7 = (red curve) and 7 = 0.25 (black curve). 
For the typical mixing strength we observed (j3 < 0.5), 
C(P) is almost the same in both cases. We can also 
notice that C(/3) increases linearly until P =0.85 and 
then decreases very fast. This is due to the fact that 
when P tends to 1, we have pure LH states and the 
total emission intensity (1% +/f) is then linearly polar- 
ized. The maximum of C(/3) does not occur for j3 = 0.5 
because as we can see in appendix A, the expressions 
of u v ±3 ^ 2 and u v ±x j 2 differ by a normalization constant 

(1/V2 and l/y/6 respectively). On the contrary, the 



Figure 7: (a) Calculated in-plane polar diagrams under non 
resonant excitation with /3 = 0.25, \j/ = 25° and 7=0 (red, 
blue and purple curves) and with 7 = 0.25 (dashed black 
curve), (b) polarization rate C(/3 , 7) versus the mixing param- 
eter p when 7 is fixed: 7=0 (dashed red curve) and 7= 0.25 
(dashed black curve). 




(a) 7 < 0.05 (b) 7=0.25 



Figure 8: Experimental polar diagrams of two dots located in 
the same detection spot, under non resonant excitation, in the 
(x,y) plane (purple) and in the (y,z) plane (orange), fitted by 
equations (6) and (12). (a) y/ = -25°, J3 = 0.2 and 7 < 0.05; 
(b) yr = -25°, P = 0.09 and 7= 0.25. 

polarization degree in the (y,z) plane greatly depends 
on 7. The transitions dipole matrix elements taking into 
account the new VBM parameter 7, give the PL nor- 
malized intensity detected along the WG x-axis I r nr (oc), 
which reads: 

4(«) = [(l--f- 

/l _ g2_y2 

+2)8 y ^ — — cos(2y)) sin 2 a 

+y 2 (- sin 2 a + ^ cos 2 a)] (12) 

The mixing parameter 7 can be determined only if 
in-plane polarization resolved jlPL experiments have 
been previously performed on the same QD to deter- 
mine P and y/. These combined experiments are diffi- 
cult to perform because we have to observe the emis- 
sion of the same QD detected from the top surface and 
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from the edge of the WG. In most cases, y has a value 
too small (y < 0.05) to be accurately estimated (Fig. 
8(a)), but for some dots 7 values range between 0.1, 
and 0.25 (Fig. 8(b)). Since 7 is related to the shear 
strain components % and £ yz , we can conclude that 
they modify also the dots optical properties and are 
partly responsible for the presence of emission polar- 
ized along the growth axis of the dots. 



V. CONCLUSION 

In summary, the polarization-resolved flPL experi- 
ments performed on single self-assembled InAs/GaAs 
QDs have allowed to determine the influence of the 
VBM on the polarization properties of the exciton 
eigenstates. Besides the strong polarization in-plane 
anisotropy, the most striking feature is that the lin- 
ear states are neither aligned along the crystallographic 
axes nor orthogonal between them. The VBM strength 
can vary significantly depending on the dot position 
along the waveguide axis while the emission polariza- 
tion axis has a given value within the detection spot 
size of one micron. Although strain and confinement 
anisotropy contributions to the VBM are difficult to 
separate, we can suggest that the strength of the mixing 
depends more strongly on confinement while the main 
polarization axis is mostly related to strain anisotropy. 
Combining the in-plane polarization-resolved jlPL ex- 
periments with resonant and off-resonance excitation, 
the ratio of oscillator strengths has been determined. 
Knowing all the polarization properties and oscillator 
strength of the linear exciton states it is then possible 
to perform real coherent manipulation of single qu-bits 
and consider future possibilities for their entanglement. 
The waveguiding geometry has also allowed observing 
emission polarized along the growth axis which is not 
possible to detect in usual backscattering photolumi- 
nescence experiments. This property is clearly related 
to shear components of the deformation tensor out of 
the growth plane. We did not find any relationship be- 
tween the strength of the mixing or the tilt of the eigen- 
states polarization axis and the emission energy of the 
exciton. We believe that the correlation between emis- 
sion energy and size of the dot can not be predicted in 
a straightforward manner especially when studying a 
large statistical ensemble of individual dots. Therefore 
it is difficult to establish general trends for the exci- 
ton fine structure versus the dots size unless an exter- 
nal field, like an electric field, 27 or stress applied to the 
structure, 17 ' 28 is used to tune the fine structure splitting. 



Appendix A: Wavef unctions 

We first assume that the band-gap energy, as well as 
the spin-orbital split-off energy are large enough to pre- 



vent any sizable coupling between the 1^, Tj and Tg 
bands, and we consider that the initial and final elec- 
tron states are eigenstates separately described by the 
valence-band Tg and the conduction-band effective- 
mass Hamiltonians 26 . In the envelope-function ap- 
proximation, the conduction band and the valence- 
band T 8 wave functions are written as: 

V /C ±i/ 2 ( r )=Xe(r)u c ±l/2 (r) (Al) 
= Xhh(r)u v ±3/2 (r) (A2) 

^ 1/2 (r)=^(r)/ ±1/2 (r) (A3) 

where % e , %hh an d Xlh are the envelope-functions and 
M ±i/2' w ±i/2' w ±i/2 are ^ e ^6 anc * ^8 Bloch functions, 
defined as 26 : 

u c 1/2 = i\st) (A4) 
u c _ x i2 = i\s \) (A5) 

«5 /2 = ^|(*+/y)t> (A6) 
K-1/2 = - ^ I (* - iy) t) - \J\ \z I) (A7) 

"V2 = ^l(*+'>H>-yf kt> (A8) 

ul 3/2 = -^=\(x-iy)i) (A9) 

where \s), \x), \y) and \z) are the orbital-functions with 
symmetry s, x, y and z, and are the spin com- 
ponents, quantized along the z-axis. For simplicity, 
we will assume that the ratio between the overlap in- 
tegrals of the electron-HH and electron-LH envelope- 
functions is equal, i.e. (Xhh\Xe) / (Xlh\Xe) ~ 1. 

Appendix B: Hamiltionian with Valence band mixing 

We assume that the valence band states in a strained 
InAs/GaAs QD can be described by a (4x4) Hamilto- 
nian based on Luttinger-Kohn and Bir-Pikus theories 
which reads 19 : 

fP + Q + V R -VlS \ 

R 1 P-Q + V V2S 
-V2S^ P-Q + V R 

\ y/2rf P + Q + Vj 

(Bl) 
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in the { Uy 2 \ tf_ij 2 , u \/2> w -3/2 } ^ as i s - 

V represents the confinement potential that usually 
can be partially separated because the confinement 
along the growth z-axis is more important than the in- 
plane confinement: V = V z (z) +V xy (x,y). The other 
matrix elements are given as a sum of kinetic terms 
and strain counterparts: 



2m 

~\~dy{£xx H~ £yy 



71 + A?) 



(B2) 



Q 



2mo 



72(k 2 x 



(B3) 



R 



+- 



V3 



\f3{Y 2 {k 2 x -k 2 y )-2ij i k x k y ) 
b v (e xx -£yy)-id v e xy (B4) 



- iky)k z 



(B5) 



7i, 72 and 73 are the modified Luttinger parameters, mo 
is the free electron mass, a v is the valence-band hy- 
drostatic deformation potential, b v and d v are the shear 
deformation potentials along the [001] and [111] axes, 
and £ij are the components of the deformation tensor. 

Using the HH band as the origin of the energies in 
the valence-band, the above Hamiltonian can be rewrit- 
ten according to Ref. 13 as: 



/ 

p se 2m 

O s 





A/a 




2*Vo 



O s 



A/a 







\ 



p s e 



-2i\j/ 



(B6) 



where A/a is the HH-LH energy separation, p s is the 
in-plane (x,y) HH-LH mixing strength, <J S is the (y,z) 
plane HH-LH mixing strength, and y/ the angle of 
in-plane anisotropy main axis with respect to [1 10]. 



Appendix C: Optical transitions matrix elements 

Without VBM, the dipole matrix element of a transi- 
tion between the conduction-band and the valence- 
band Tg states can be written as: 

<V^|e.p|v4> = <±l|e-P|0> 

= (Xe\Xhh) (u C Tl/2 \£.p\u V ±3/2 ) (CI) 

where the polarization vector £ defines the orientation 
of the electric field of the linearly polarized excitation 
laser and p is the electron linear momentum. When 
VBM exists, the valence-band states become: 



m) = Vi-P 2 -r 2 \x hh )w ± y 2 ) 

= +/^ y l*»>Ki/2> 

+ye ±2 * \Xih) K ±1/2 ) (C2) 

and the dipole transition matrix element between con- 
duction and valence states is modified and reads: 

<V^|e.p|^> = (±l|e.p|0> 

= (Xe\Xhh) V^-P 2 -f(u c Tl/2 \£^\u v ±3/2 ) 

+ (^|to)^^ >0 «l/2l^Pl<l/2) 

+ (Xe\Xih)y(u c T y 2 \£.v\u v ± 1/2 ) (C3) 

We assume that the overlaps (Xhh\Xe) / (Xih\Xe) ~ 1, 
and we introduce the mixing parameters j3 and 7, re- 
lated to p s , o s and A ih by: A ^ 9 = = ^ and 
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